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Abstract 

Determinations of the strong coupling strength, a,, are summarised and a new world average value of OsiMz) is 
determined, using a new method of pre-averaging results within classes of measurements like hadronic r decays, deep 
inelastic scattering processes, lattice calculations, electron-positron annihilation processes and electro- weak precision 
fits. The overall result is 

ff,(Mz) = 0.1 184 ±0.0007, 

unchanged from the value obtained in 2009. This presentation is an excerpt from the QCD review section of the 2012 
Review of Particle Physics (RPP) of the Particle Data Group yj. An earlier version of this work was also given in 10]. 
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1. Introduction 

The size of the strong coupling strength, as, like other 
fundamental "constants" of nature, is not given by the- 
ory, but must be determined by experiment. Quantum 
Chromodynamics (QCD) the gauge field theory de- 
scribing the Strong Interaction between quark and glu- 
ons, the fundamental constituents of hadronic matter, al- 
lows to predict physical cross sections and many other 
observables K of particle reactions involving quarks and 
gluons, in terms of - basically - one single parameter, 
as. Assuming as < 1 and applying perturbation theory, 
predictions are typically given by a power series in as, 
hke: 



= P,[R(,+R,a,+R2al + ..) (1) 

where R„ are the n,h order coefficients of the perturba- 
tion series and PiRq denotes the lowest-order value of H. 
This allows to measure as in a large variety of particle 
reactions. 

QCD also predicts the energy dependence of as 
through the "beta function". 



(2) 
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where is the squared energy scale or the momen- 
tum transferred of the particle reaction under study, and 
any energy scale dependence of H is determined by the 
energy dependence of afl Measurements of as from 
different particle reactions and scattering processes, per- 
formed at different energy scales Q^, therefore test the 
global nature of QCD and, in particular, its character- 
istic prediction of "Asymptotic Freedom", which deter- 
mines that as is small and asymptotically tends to zero 
at large energy scales (or at small distances), while it is 
large at small energy scales (large distances), explaining 
the "confinement" of quarks and gluons inside hadrons. 

Experimental determinations of as were regularly 
summarised and reviewed in the past, see e.g. 
These references also contain the definition of basic 
equations and fomulae which shall not be repeated here. 
In 2009, the world average value of as, expressed at a 
common energy scale corresponding to the rest mass of 
the Z boson, determined from a set of most recent de- 
terminations from many different processes at a large 
range of energy scales, converged to 

as(Mz) = 0.1184 + 0.0007, 

where the overall error includes experimental as well 
as (dominating) systematic and theoretical uncertainies 



' For processes with initial state hadrons, also the coefficients i?„ 
may depend on a,, through the effects of parton density functions. 



Preprint submitted to Niic. Phys. (Proc. Suppl.) 



October 2, 2012 



In the following, this summary is updated with the 
newest results in this field. 

2. Selected results 

The set of results used in this review is restricted to 
those which were published in peer-reviewed journals 
until April 2012, and which use QCD perturbation the- 
ory to at least full next-to-next-to-leading order pertur- 
batioiJl . 

These requirements exclude e.g. results from jet 
production in deep inelastic lepton-nucleon scattering 
(DIS) and from hadron collisions, as well as those 
from heavy quarkonia decays for which calculations are 
available in NLO only. While being excluded from cal- 
culating a world average value of ffj, NLO results will 
nevertheless be cited in this review as they are impor- 
tant ingredients for demonstrating the experimental ev- 
idence of the energy dependence of c,, i.e. for Asymp- 
totic Freedom, one of the key features of QCD. 

Furthermore, here we add an intermediate step of 
pre-averaging results within well defined sub-fields like 
e^e~-annihilation, DIS and hadronic r-decays. The 
overall world average is then calculated from those 
pre-averages rather than from individual measurements. 
The intermediate step is done because in a number of 
sub-fields, different determinations of the strong cou- 
pling from substantially similar datasets lead to values 
of Us that are only marginally compatible with each 
other, or with the final world average value, which may 
be a reflection of the challenges of evaluating system- 
atic uncertainties. In such cases, a pre-average value 
is determined, with a symmetric, overall error that en- 
compasses the central values of all individual determi- 
nations. 

2.1. Hadronic T decays 

Determination of cKs from hadronic t lepton decays 
continues to be one of the most actively studied fields to 
measure this basic quantity. The small effective energy 
scale, Q = - 1.78 GeV, small nonperturbative con- 
tributions to an inclusive and well defined experimental 
observable, and the availability of perturbative predic- 
tions which are complete to N^LO determine the impor- 
tance and large interest in this particular field. Several 
re-analyses of the hadronic t decay width ||7|-Il2| were 
performed, using different approaches of treating per- 
turbative (fixed order or contour improved perturbative 



expansions) and non-perturbative contributions. The re- 
sult of [f] includes both, fixed order and contour im- 
proved perturbation, while the others adhere to either 
one or the other of the two. These results are sum- 
marised in Fig. [Ha). 

There are more studies of ffj from r-decays, lfl3l - [l6ll . 
which were not yet available as peer-reviewd publica- 
tions. They are compatible with the overall picture as 
summarized in Fig. [TJa). Another study |17] argues 
that an improved treatment of non-perturbative effects 
results in values of a,, which are systematically lower 
than those discussed above. These results have recently 
been extended in BlSll . 

The pre-average result from r-decays, to be used for 
calculating the final world average of q's(M|), is deter- 
mined using the simple method mentioned above, i.e. 
defining one central value with symmetric overall error 
bars which include the smallest as well as the largest 
of all results, as a,(M?) = 0.330 ± 0.014. This value 
of as{M^) corresponds, when evolved to the scale of 
the Z-boson, using the QCD 4-loop beta-function plus 
3-loop matching at the charm- and the bottom-quark 



NNLO; for observables with QCD contributions starting at lead- 
ing order, that is 0(al). 



masses (see Il4|-|6t] for relevant equations), to as(M^) = 
0.1 197 ± 0.0016, unchanged from its value in the 2009 
review. 



2.2. Lattice QCD 

There are several recent results on ff, from lattice 
QCD, see also Sec. Lattice QCD in UJ. The HPQCD 
collaboration |19] computes short distance quantities 
like Wilson loops with lattice QCD and analyzes them 
with NNLO perturbative QCD. This yields a value for 
a,. The lattice scale must then be related to a phys- 
ical energy/momentum scale. This is achieved with 
the T'-T mass difference, however, many other quan- 
tities could be used as well 1I20I1 . HPQCD obtains 
a, (Ml) - 0.1184 ± 0.0006, where the uncertainty in- 
cludes effects from truncating perturbation theory, finite 
lattice spacing and extrapolation of lattice data. An in- 
dependent perturbative analysis of a subset of the same 
lattice-QCD data yields a',(M|) = 0.1192±0.0011||21|. 
Using another, independent methodology, the current- 
current correlator method, HPQCD obtains q's(M|) = 
0.1183 ± 0.0007 [19]. A more recent result in 121, 
which avoids the staggered fermion treatment of |19], 
finds a,(Ml) = 0.1205 ± 0.0008 ± 0.0005 !{|j;°°°r2l, 
where the first uncertainty is statistical and the others 
are from systematics. Since this approach uses a dif- 
ferent discretization of lattice fermions and a different 
general methodology, it provides an independent cross 
check of other lattice extractions of a,. Finally, the 
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JLQCD collaboration - in an analysis of Adler functions 
- obtains a,(Ml) = 0.1 181 + 0.0003 !„:„oi'2l'23']. 

These results are summarized in Fig. [Tib). Since they 
are compatible with and largely independent from each 
other, a pre-average of lattice results is calculated us- 
ing the same method as applied to determine the final 
world average value as, i.e. calculate a weighted av- 
erage and a (correlated) error such that the overall 
equals unity per degree of freedom - rather than using 
the simple method as applied in the case of r decays. 
This gives a,(Mz) = 0.1185 ± 0.0007 which is taken as 
result from the sub-field of lattice determinations. 



2.3. Deep inelastic lepton-nucleon scattering (DIS) 

Studies of DIS final states have led to a number of 
precise determinations of as: 

A combination |24] of precision measurements at 
HERA, based on NLO fits to inclusive jet cross sections 
in neutral current DIS at high Q^, quotes a combined 
result of a,(Mz) = 0.1198 + 0.0032, which includes a 
theoretical uncertainty of ±0.0026. A combined anal- 
ysis of non-singlet structure functions from DIS ll25ll . 
based on QCD predictions up to N^'LO in some of its 
parts, gave a^iMz) - 0.1 142 ± 0.0023, including a the- 
oretical error of ±0.0008. Further studies of singlet and 
non-singlet structure functions, based on NNLO predic- 
tions, resulted in a,(Mz) = 0.1129 ± 0.0014 [26] and 
in as(Mz) = 0.1158 ± 0.0035 t27J. The MSTW group 
0281], also including data on jet production at the Teva- 
tron, obtains, in NNLcH a,(Mz) = 0.1171 ± 0.0014. 
Most recently, the NNPDF group 1E9I1 has presented a 
result, as(Mz) = 0.1 173 ± 0.001 1, which is in Une with 
the one from the MSTW group. 

Summarizing these results from world data on struc- 
ture functions, applying the same method as in the 
case of summarizing results from r decays, leads to a 
pre-average value of a^iMz) - 0.1151 ± 0.0022 (see 
Fig.Qc)). 

Note that criticism has been expressed on some of the 
above extractions. Among the issues raised, we mention 
the neglect of singlet contributions at x > 0.3 in pure 
non-singlet fits 13011 . the impact and detailed treatment 
of particular classes of data in the fits and possi- 

ble biases due to insufficiently flexible parametrizations 
ofthePDFs Q. 



^^Note that for jet production at tlie hadron collider, only NLO pre- 
dictions are available, while for the structure functions full NNLO was 
utilized. 



2.4. Heavy quarkonia decays 

The most recent extraction of the strong coupling 
constant from an analysis of radiative T decays [33'] re- 
sulted in a^iMz) - 0.1 19+J5 "°g. This determination is 
based on QCD in NLO only, so it will not be considered 
for the final extraction of the world average value of a/, 
it is, however, an important ingredient for the demon- 
stration of Asymptotic Freedom as given in Fig. [3] 

2.5. Hadronic final states ofe^e^ annihilations 
Re-analyses of event shapes in e^e -annihilation, 

measured at the Z peak and LEP2 energies up to 209 
GeV, using NNLO predictions matched to NLL resum- 
mation, resulted in aJMz) = 0.1224 ± 0.0039 [3, 
with a dominant theoretical uncertainty of 0.0035, and 
in as(Mz) = 0.1189 + 0.0043 Q. Similarly, an anal- 
ysis of JADE data f36'l at center-of-mass energies be- 
tween 14 and 46 GeV gives a^iMz) = 0.1 172 ± 0.0051, 
with contributions from hadronization model (pertur- 
bative QCD) uncertainties of 0.0035 (0.0030). A pre- 
cise determination of from 3-jet production alone, 
in NNLO, resulted in aJMz) = 0.1175 ± 0.0025 Q. 
Computation of the NLO corrections to 5 -jet production 
and comparison to the measured 5-jet rates at LEP Ii 38ll 
gave as(Mz) - O.II56+QQ034. More recently, a study 
using the world data of Thrust distributions and soft- 
coUinear effective theory, including fixed order NNLO, 
gave a,(Mz) = 0.1135 ± 0.0010 |39l]- 

Note that there is criticism on both classes of a, 
extractions just described: those based on corrections 
of non-perturbative hadronisation effects using QCD- 
inspired Monte Carlo generators (since the parton level 
of a Monte Carlo is not defined in a manner equivalent 
to that of a fixed-order calculation), as well as the stud- 
ies based on effective field theory, as their systematics 
have not yet been verified e.g. by using observables 
other than Thrust. 

A summary of the e'^e' results based on NNLO pre- 
dictions is shown in Fig. [Hd). They average, accord- 
ing to the simple procedure defined above, to ci's(Mz) - 
0.1172 ±0.0037. 

2.6. Hadron collider jets 

A determination of a, from the pr dependence of 
the inclusive jet cross section in pp collisions at = 
1 .96 TeV, in the transverse momentum range of 50 < 
Pt < 145 GeV, based on NLO (Oiaj)) QCD, led to 
UsiMz) = 0.1161+[|;[|°^i [40J, which is the most precise 
as result obtained at a hadron collider. Experimental un- 
certainties from the jet energy calibration, the pr resolu- 
tion and the integrated luminosity dominate the overall 
error. 
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Figure 1: Summary of determinations of a, from hadronic r-decays (a), from lattice calculations (b), from DIS structure functions (c) and from 
event sliapes and jet production in e+e^-annihilation (d). Tlie sliaded bands indicate the average values chosen to be included in the determination 
of the new world average of qtj. Figure taken from [IJ. 



2.7. Electroweak precision fits 

The N^LO calculation of the hadronic Z decay width 
was used in a recent revision of the global fit to elec- 
troweak precision data 114 ill , resulting in as{Mz) - 
0.1193 + 0.0028, claiming a negligible theoretical un- 
certainty. For this review, we use the result obtained in 
H from data at the Z-pole, a,(Ml) = 0.1197 ± 0.0028, 
as it is based on a more constrained data set where QCD 
corrections directly enter through the hadronic decay 
width of the Z. Note that all these results from elec- 
troweak precision data, however, strongly depend on the 
strict validity of Standard Model predictions and the ex- 
istence of the minimal Higgs mechanism to implement 
electroweak symmetry breaking. Any - even small - de- 
viation of nature from this model could strongly influ- 
ence this extraction of a^. 

3. Determination of the world average value of 

A non-trivial exercise consists in the evaluation of a 
world-average value for a^(Mz)- A certain arbitrariness 



and subjective component is inevitable because of the 
choice of measurements to be included in the average, 
the treatment of (non-Gaussian) systematic uncertain- 
ties of mostly theoretical nature, as well as the treat- 
ment of correlations among the various inputs, of theo- 
retical as well as experimental origin. In earlier reviews 
13-0 an attempt was made to take account of such cor- 
relations, using methods as proposed e.g. in (a^], and 
- likewise - to treat cases of apparent incompatibilities 
or possibly underestimated systematic uncertainties in a 
meaningful and well defined manner: 

The central value is determined as the weighted av- 
erage of the different input values. An initial error of 
the central value is determined treating the uncertain- 
ties of all individual measurements as being uncorre- 
lated and being of Gaussian nature, and the overall 
to the central value is determined. In case this initial 

is larger than the number of degrees of freedom, i.e. 
larger than the number of individual inputs minus one, 
then all individual errors are enlarged by a common fac- 
tor such that /d-o.f. equals unity. If the initial value 
ofx^ is smaller than the number of degrees of freedom. 
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an overall, a-priori unknown correlation coefficient is 
introduced and determined by requiring that the total 
Id.o.f. of the combination equals unity. In both cases, 
the resulting final overall uncertainty of the central value 
of Us is larger than the initial estimate of a Gaussian er- 
ror. 

This procedure is only meaningful if the individual 
measurements are known not to be correlated to large 
degrees, i.e. if they are not - for instance - based on 
the same input data, and if the input values are largely 
compatible with each other and with the resulting cen- 
tral value, within their assigned uncertainties. The list 
of selected individual measurements discussed above, 
however, violates both these requirements: there are 
several measurements based on (partly or fully) iden- 
tical data sets, and there are results which apparently do 
not agree with others and/or with the resulting central 
value, within their assigned individual uncertainty. Ex- 
amples for the first case are results from the hadronic 
width of the t lepton, from DIS processes and from jets 
and event shapes in e^e" final states. An example of the 
second case is the apparent disagreement between re- 
sults from the t width and those from DIS fzS*] or from 
Thrust distributions in e^e" annihilation [39]. 

Due to these obstacles, we have chosen to determine 
pre-averages for each class of measurements, and then 
to combine those to the final world average value of 
OsiMx), using the methods of error treatment as just de- 
scribed. The five pre-averages are summarized in Fig.|2j 
we recall that these are exclusively obtained from ex- 
tractions which are based on (at least) full NNLO QCD 
predictions, and are published in peer-reviewed journals 
at the time of completing this review. From these, the 
new central and world average value of 

a,(Mz) = 0.1184 ±0.0007, (3) 

is determined, with an uncertainty of well below 1 %0 
This world average value is - in spite of several new 
contributions to this determination - identical to and 
thus, in excellent agreement with the 2009 result 
For convenience, we also provide corresponding val- 
ues for h-f^ suitable for use with the common A- 
parametrisation of a^, see e.g. Eq. 6 in [5]: 

Ai!L = (213±8)MeV, (4) 

MS ^ ' 

AilL = (296+ 10)MeV, (5) 



"^The weighted average, treating all inputs as uncorrelated mea- 
surements with Gaussian errors, results in a^{Mz) = 0.11844 ± 
0.00059 with ;('^/d.o.f. = 3.2/4. Requiring ^^/d.o.f. to reach unity 
leads to a common correlation factor of 0. 19 which increases the over- 
all error to 0.00072. 
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Figure 2: Summary of values of a^iMz)) obtained for various sub- 
classes of measurements (see Fig.[T](a) to (d)). The new central and 
world average value oi a^{Mz) = 0.1184 ± 0.0007 is indicated by the 
dashed line and the shaded band. Figure taken from \ \^. 

h2L = (339 ± 10) MeV , (6) 

MS ^ ' 

for Nf = 5,4 and 3 quark flavors, respectively. 

In order to verify the consistency and stability of the 
new result, we give each of the averages obtained when 
leaving out one of the five input values: 

a,(Mz) = 0.1182 + 0.0007 (w/o r results), 

a,(Mz) = 0.1183 + 0.0012 (w/o lattice), 

q;,(Mz) = 0.1187 + 0.0009 (w/o DIS), 

as(Mz) = 0.1185 + 0.0006 (w/o e+el, and 

as(Mz) = 0.1184 + 0.0006 (w/o e.w.prec. fit). 

They are well within the error of the overall world av- 
erage quoted above. Most notably, the result from lat- 
tice calculations, which has the smallest assigned error, 
agrees well with the exclusive average of the other re- 
sults. However, it largely determines the size of the 
(small) overall uncertainty. 

There are apparent systematic differences between 
the various structure function results, and also between 
the new result from Thrust in e^e" annihilation and the 
other determinations. Expressing this in terms of a 
between a given measurement and the world average as 
obtained when excluding that particular measurement, 
the largest values are x^ - 12.6 and x^ - 16.1, cor- 
responding to 3.5 and 4.0 standard deviations, for the 
measurements of [2^ and ^39'], respectively. We note 
that such and other differences have been extensively 
discussed at a recent workshop on measurements of a,,, 
however none of the explanations proposed so far have 
obtained enough of a consensus to definitely resolve 
these tensions 04311 . 
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Figure 3: Summary of measurements of a, as a function of the re- 
spective energy scale Q. The respective degree of QCD perturbation 
theory used in the extraction of ff , is indicated in brackets (NLO: next- 
to-leading order; NNLO: next-to-next-to leading order; res. NNLO: 
NNLO matched with resummed next-to-leading logs; N'LO: next-to- 
NNLO). Figure taken from 

Notwithstanding these open issues, a rather stable 
and well defined world average value emerges from the 
compilation of current determinations of a,: 

a,(Mz) = 0.1184 + 0.0007 . 

The results also provide a clear signature and proof of 
the energy dependence of a,, in full agreement with 
the QCD prediction of Asymptotic Freedom. This is 
demonstrated in Fig. [3] where results of OsiQ^) obtained 
at discrete energy scales Q, now also including those 
based just on NLO QCD, are summarized. 
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